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We show that under particular circumstances a general relativistic spherically symmetric bounded 
distribution of matter could satisfy a nonlocal equation of state. This equation relates, at a given 
point, the components of the corresponding energy momentum tensor not only as function at that 
point, but as a functional throughout the enclosed configuration. We have found that these types 
of dynamic bounded matter configurations, with constant gravitational potentials at the surface, 
admit a Conformal Killing Vector and fulfill the energy conditions for anisotropic imperfect fluids. 
We also present several analytical and numerical models satisfying these equations of state which 
collapse as reasonable radiating anisotropic spheres in general relativity. 



I. INTRODUCTION 



Our comprehension of the behavior of highly compact stars is intimately related to the understanding of the physics 
at supranuclear densities. Today, the properties of matter at densities higher than nuclear ( « lO^** gr/cm.^) are 
essentially unknown, although they must be oriented by the experimental insight and extrapolations emerging from 
the ultra high energy accelerators and firmly anchored to what is known from nuclear physics Having this 

uncertainty in mind, we shall explore what is allowed by the laws of physics, for a particular equation of state, within 
the framework of the theory of General Relativity and considering spherical symmetry. 

We shall consider a spherically symmetric space-time that can be described by the metric 

ds^ = h e^^dT^ - \dR^ - R^dn^ , (1) 
h 
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where the sohd angle is dfi^ = dO^ + sin^ 6d(jy^ and both /3 and h depend on T and R. 
We have found that if in the above space-time (|l|) an additional restriction, 

/i(r,i?) = i-^!^%^ = c(T)e-2/?(T,i?)^ ^-^^ o<c(r)<i, (2) 

is considered, we can obtain the following relation between two of the components of the corresponding Einstein 
tensor, i.e. 

G? + 3Gg + i?A[G? + Gg] =0. (3) 
Equation (||) can be integrated yielding 



G^ = -Gj^--y^Gi^ df + ^ , (4) 
which, using the Einstein equations, can be re-written as 

where 1C{T) is an arbitrary integration function. Obviously, from (^, we can also solve Gj^ in order to obtain an 
expression equivalent to (ra), i.e. 



Again, ^(T) is an arbitrary integration function. Clearly, for bounded matter distributions both functions K,{T) and 
G{T), should emerge from the junction conditions to match the corresponding exterior metric. 

It can be easily appreciated that the above equations (|^) and involve a non local relation between the components 
of the energy momentum tensor. At a given instant of time, one of the component (say in (|^)) is not only a 
function of the other ( ) at that point but also its functional throughout the rest of the configuration. In the case 
of radiating fluids, equations (||) and (||) relate non-locally the hydrodynamic and radiation physical variables. 

Equation (^) could be further re-written as 

T^ = -T? + ^(T?)-^, (7) 

with 

(^^t) - Y¥ jJ^^r^T^ df , (8) 

clearly the nonlocal term represents an average of the function T'^^ over the volume enclosed by the radius R. Moreover, 
equation can be easily rearranged as 

where we have used the concept of statistical standard deviation ar^r from the local value of the function TJ! . 
Therefore, if at a particular point within the distribution the value of the function gets very close to its average 
(Tx) relation between and T'^ becomes 

« -^T? - ^ with « . (10) 

Physical insight can be gained by considering equations (^ and (^) in their static limit and in terms of the 
hydrodynamic density, p{R) and pressure, Pr{R). In this case, a static spherically symmetric distribution leads to 
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p(i?) = P,(i?) + |y^P,(f) df + ^ , or (11) 



Pr{R) = p[R) - ^ f'p{f) df + ^ 



(12) 



with H and B are arbitrary integration constants. It is clear that in equations (|Tl|) and (|T^ a coUective behavior on 
the physical variables p{R) and Pr{R) is also present, and could be interpreted as a nonlocal interrelation between the 
energy density, p{R), and the hydrodynamic pressure, Pr{R), within the fluid. Any change in the pressure (density) 
takes into account the effects of the variations of the energy density (pressure) within an entire volume. 
Following the above line of reasoning we shall consider the static limit of the equation (^: 

P,(i?)=7'(i?) + 2ap(^) + -^ , (13) 

where we have set 

V{R) = l^piR) and a-piR) = Qp(i?) - ^ (p)^) = {V{R) - V{R)) , (14) 



(P). = ^f 4..-.(^-)d.=|f. (15) 



with 



In this limit the nonlocal integral term in ( p_2D represents some kind of "average density" over the enclosed volume 
and, as far as the value of this nonlocal contribution gets closer to the value of the local density, the equation of state 
of the material becomes similar to the typical radiation dominated environment, i.e. Pr{R) ~ "PiR) = ^p{R)- 

It is in the above sense, that we are going to refer to (^ or to (^, as Nonlocal Equation of State {NLES from now 
on) between the components and T|| of the corresponding energy momentum tensor. 

Nowadays, the relevance of nonlocal outcomes on the mechanical properties of fluids and materials are well known. 
There are many situations of common occurrence wherein nonlocal effects dominate the macroscopic behavior of 
matter in Modern Classical Continuum Mechanics and Fluid Dynamics. Examples coming from a wide variety of 
classical areas such as: damage and cracking analysis of materials, surface phenomena between two liquids or two 
phases, mechanics of liquid crystals, blood flow, dynamics of colloidal suspensions seem to demand more sophisticated 
continuum theories which can take into account different nonlocal effects. This is a very active area in recent material 
and fluid science and engineering (see |^ and references therein). Moreover, in radiating fluids when the medium 
is very transparent and photons escape efficiently, occupation numbers of atomic levels are, in general, no longer 
predicted by equilibrium statistical mechanics for local values of the temperature and density. In this case, the local 
state of a fluids is coupled by photon exchange to the state of the material within an entire interaction volume. These 
types of fluids are considered to be in non-Local Thermodynamical Equilibrium or non-LTE Q . 

In this work only collapsing radiating configuration in General Relativity will be studied. Static models consistent 
with NLES require a detail analysis and will be considered elsewhere ||]. In what follows, we shall explore four major 
questions concerning the NLES in collapsing spherical sources in General Relativity, which can be stated as follows: 

• Can matter configurations having NLES be associated with a particular symmetry of the space time and under 
what conditions ? 



Does the energy- momentum tensor obtained from (^ through the Einstein Equations, and satisfying (^, describe 
a physically reasonable anisotropic fluid ? 

• Can this fluid represent a bounded gravitational source ? 

• What types of gravitational collapse scenarios will emerge from bounded matter configurations with NLES 7 

The structure of the present paper is guided by the above four questions. In the next section it is found that 
a solution of the Einstein equations, obtained from (^ restricted by (H), admits a Conformal Killing Vector. In 
Section III it is shown how anisotropic radiating configurations with NLES satisfy energy conditions (weak, strong 
and dominant) for imperfect fluids. The third question is addressed in Section IV. There we consider the consequences 
on the evolution of the boundary surface when this interior solution is matched to Vaidya exterior metric. Answers 
to the last question are explored in sections IV and V by studying the evolution of models of radiating spheres with 
a nonlocal equation of state (either |^ or Finally in the last section our conclusions and results are summarized. If 
it is not explicitly stated as C(T) (or C(u), depending on the tetrad) we shall assume in this paper C as a constant 
parameter in the above equation (g). 
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II. SYMMETRIES AND NLES FOR AN ANISOTROPIC RADIATING SPHERE 



It is well known that there are, at least, two reasons why it is interesting to start investigating the geometric 
properties of a particular space time. The first is a practical one: the formidable difhculties encountered in the study 
of the Einstein Field Equations, due to their non-linearity, can be partially overcame if some symmetry properties of 
the space-time are assumed. The second one is more fundamental; symmetries may also provide important insight and 
information into the general properties of self-gravitating matter configurations [|6| . Guided by the above motivations 
we shall show, in this section, that the space times having a NLES or|6|) admit a Conformal Motion [CM) with 
C considered as a constant parameter. 

In general, CM is a map M ^ M such that the metric g of the space time transforms under the rule 



with 



This can be expressed as 



£i gab = £.a;b + ib;a = i^i^^) gafc , 



(16) 



where is called a Conformal Killing Vector [CKV). It is clear that other important symmetries as homotetic motions 
or self-similarity gab — dab with a — const.) and isometrics gab = 0) are particular cases of CM. This 
symmetry imposes an important restriction on the hydrodynamic variables and consequently it is possible to obtain 
an equation of state for a space time (|l|) having a CKV [7|-pX|].For the sake of simplicity and in order to solve the above 



Conformal Killing Equations (|16|) we shall assume C as a constant parameter. The most general case, C = C(T), will 
be briefly considered at the end of the present work. 

Additionally, we further restrict our calculations to a vector field of the form 



From m and <\l 



r = CT(T,i?K + A(T,flX 

([l^), with the condition (||), we have 



A' 



A -f cr/3 = , 



-A + C^ct' = 0. 



and 



A -f cr/3 = . 



(17) 

(18) 
(19) 
(20) 



where dots and primes denote differentiation with respect to T and i?, respectively. 
Using ^ and (|^) we get 

-A + C2A" = 0. (21) 
A similar expression can be obtained for a. The solutions of equation ( pl|) and the corresponding equation for cr are 

A = /(u) -f g(v) , and cr = m(u) + n{v) , (22) 
/(u), g(y), m(u) and »t.(v) being arbitrary functions of their arguments 

u = CT -I- i? , and v = CT - i? . (23) 
With equations (|||) in (|l|), we obtain 

1 



C 



[/(n)-5(v)] 



and, form (G2h and 



equation (20) can be written as 
1 



d(3_ 



dn 



dp 
dv 



dg „ 



(24) 



(25) 



This can be integrated by demanding that the coefficient of /(u) be only a function of u as well as the coefficient g{v) 
a function of v: 
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da{u) _ d(3 



and 



dS{v) _ dp 



1 



dn 9u u — V ' dv (9v u — v 

Thus, we get 

/3(u, v) = In (u - v) + q;(u) + 5{v) + Ci , 
and equation ( p5| ) may now be written as 

„ fda\ , df fdS\ dg , 

where k and Ci are constants of integration. This expression is integrated to give 



/(u) = e- 



k I e2"(iu+ C2 



and g{v) = e 
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k I e^^dY+C3 



again, C2 and C3 are constants of integration. Then, ( p9| ) impHes 

e = ^ [/(u) - g(v)] + [/(u) + g(v)] and [/(u) + ^(v)] , 

and ( p^ ) may now be written as 

/3(T, i?) = ln(2fci?) + a(CT + i?) + (5(CT - R) + Ci . 

Finally, the metric (El) reads: 



(26) 



(27) 



(28) 



(29) 



(30) 



(31) 



(32) 



It should be stressed that the above results have been obtained assuming (E9) and C as a constant parameter in 



III. ENERGY CONDITIONS FOR IMPERFECT FLUIDS 



We have found that in static configurations NLES are more easily handled within anisotropic fluids ||]. Thus, in 
order to explore how physically reasonable is the previous non-static fluid with a NLES and having a CKV, we shall 
study the energy conditions for an imperfect anisotropic fluid (unequal stresses, i.e. Pr ^ P± and heat flux). Although 
a pascalian assumption {Pr = P±) is supported by solid observational ground, an increasing amount of theoretical 
evidence strongly suggests that, for certain density ranges, a variety of very interesting physical phenomena may take 
place giving rise to local anisotropy (see [prf and references therein). 

For these fluids the energy-momentum tensor takes the form 



T^<. = (P + P±)u^u^ - Pig^,. + {Pr - P±)n^n^ + f^Ui, + f^u^ 



where 



u 



and 



fu = -q 



(33) 



(34) 



with p, Pr, P± and q, denoting the energy density, the radial pressure, the tangential pressure and the heat flow 
(f^jf = —9^), respectively. The fluid four- velocity is represented by and is a unit vector pointing the flux 
direction. 



The energy conditions are obtained from the eigenvalues of the energy-momentum tensor (33), i.e. they emerge 
from the roots of the characteristic equation 



(35) 



The four eigenvalues take the form 



5 



Ao = ^ [p - Pr + A] , Ai = i [p - p. - A] , and A2 = A3 = -Pj. , (36) 

where 

A=[(p + P.)^-V]l/^ (37) 
With the metric (^2[), Einstein's field equations can be written as 

87rT? = 8^p=i2 +-^[l + 2P(a'-<5')] , (38) 
87rTg = -87rP, = ij--^[3+ 2R{a' - 5')] , (39) 

87rTg = -87rP_i = and (40) 

-ft 

87rT™ = -87r<7= — + , (41) 

where 



From two of the above Einstein equations, (|3^) and (p9[), we find that the conditions for energy density and radial 
pressure to be positive are 

a'-5'>-—-— and a' ^ 5' > — - — . (43) 

Considering (^0|) and ( ^2[ ) it can easily be shown that the tangential pressure, P±, is negative. Consequently the 
condition (^7|) is fulfilled if 

2P_i + A > ^ a' -S' >- (a + s) - - . (44) 



P 



We then find that the weak energy condition is achieved if: 



2Z 2 

p-Pr + A>0 ^ A>— - — , (45) 



the dominant energy condition is held if 



Z 2 

p-P^>0^— r<l and p-P^-2Pj. + A>0 ^ A> r, (46) 
P^ R-' 

and the strong energy condition is accomplished if 

2Z 

2P_i + A>0=»A> — . (47) 

Next, we are going to show that it is possible to fulfill simultaneously all of the above energy conditions. From the 
definition for Z and (^) it is now seen that 

0<J<1. (48) 

Therefore, inequalities (^), (^) and (|4^) yield 



2 - P2 - 2AP4 

Now, considering (^), (^) and (^4|) we have 

1 



> > (49) 



d + (5>0. (50) 
P 



By using equation (pn) we can also see that 



P' = ^+a'-S', (51) 



which, with the help of ( pO| ) becomes 

P'>0. (52) 

This is a well known restriction on the metric function (3 [p^ . Therefore, it is clear that matter configurations 
satisfying a NLES and having a CKV can represent reasonable fluids in General Relativity. Next section will be 
devoted to study the consequences of matching this type of fluids to an exterior solution. 



IV. JUNCTION CONDITIONS 



In order to address the third question, we shall demand that the line element ( p2| ) join the Vaidya exterior metric 
at the boundary of the configuration. We shall also recall that two regions of a space-time match across a separating 
hypersurface S if the first and the second fundamental forms are continuous across this boundary surface S |14|. 

The Vaidya line element outside the source can be written as 



dsf_^^^ = Hdu^ + 2dudr - r'^dVt^ , 



with 



H = 1 



2M (u) 



(53) 



where the subscript (+) refers to the exterior region of the space time. In these coordinates the boundary equation 
takes the form: r — rs{u), where the subscript s reminds that the corresponding function has been evaluated at 
boundary of the configuration. Then, the induced metric on the boundary surface, from the outside, is 



ds 



2\i+) 



H, 



^ drs 
du 



du^ - rUn^ . 



and the line element on the boundary, from the inside, reads 



^ 1 (dRs 



dT^ - R^dn^ , 



(54) 



(55) 



where we have used the fact that the equation of the boundary in the coordinates (T, i?, 0,<I>) reads: R = Rs{T). 
Demanding continuity of the first fundamental form, we get at once 



C 



1 / dR. 



C\dT 



It can be shown [|5| that (E^) is equivalent to 



1/2 

dT : 
Rs (T) : 



1 



) 1 

du 



1/2 



du . 



and 



(56) 
(57) 

(58) 



Also, a straighforward calculation shows that the continuity of the second fundamental form across the boundary 
surface is equivalent to the condition 



kl = [PAs ■ 

Using ( ^9|) and (|4l|), the equation (^) gives the evolution of the boundary surface: 

1 



Rs = 



-Si 



(59) 



(60) 



where the constant k is 
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K=—^. (61) 
Now, substituting equation ( |6C| ) in (|58|), we obtain 

+ + a, - 5^ = K e^^^-^^ . (62) 
It is useful to write equation (^2|) in terms of the variables u and v 

fda, + (C + 1) = «e("^(")+*'(-)) , (63) 



(64) 



\ du dv 

which can be easily integrated for constant a^. Thus, assuming = K = const, we get 

r A.(v-(c + i)C2) 

S.^-K-\n ^ 

Next, considering K and Ss in equation (|5^), the evolution of the boundary surface can be written as 

K±2{C + l) ' ^^^^ 

where C2 is a constant of integration. It is clear that there are two possible evolutions for the bounding surface and 
in both of them Rs is a linear function of the time coordinate, namely 



_ (1-3C)C _ (1-3C)(C + 1)C+ 
- 2C2 - C + 1 ^ 2C2 - C + 1 ' ^'''^ 



P-^ (1-3C)C (1-3C)(C + 1)C- 

-2C2-5C + 1 -2C2-5C + 1 ■ ^ 

Notice that expression (^l|) for the constant k, has been used and we have denoted C2 and C2 the corresponding 
constants for the two equations emerging from (|65|). 

Let us consider each of the above solutions separately. First, it is clear that in (p6|), due to 



20^ - C + 1 > V C e (0, 1) , (68) 

expanding configurations can be obtained if C e (O , |) and collapsing ones for C G ( | , l) • Additionally, because 
Rf > V T G (0,00), several values of the possible initial configurations (Mq and R^) have to be excluded. 
Concerning the first case, the allowed values for the constant Care 

C e (0, Expanding Configurations, and (69) 



C e , 1^ ^ Contracting Configurations. (70) 
In the second case (equation (|67|)), the permitted values can be expressed as 
C e (0, ifV33-5n U ( i , 1 



=^ Expanding Configurations, and (71) 
C e ^-|^V33 — 5^ , Contracting Configurations. (72) 

It is clear that C2 and C2'can be solved from (|6^), and (|67|), namely 

i?+ (2C2-C + 1) 

^ ' (1-3C)(C + 1) ' ""'^ ^''^ 
_ _ i?,- (2 C2 + 5 C - 1) 
^2 - (1-3C)(C + 1) ' ^^^^ 
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with Rfg and Rj^ the initial values for Rf and R~ , respectively. For the present work we shall only consider those 
values of C that induce collapsing configurations. 

For this particular matter distribution, starting with an initial radius, Rf^, and mass, Mq, we shall define the Total 
Evolution Time, Tev, as the time it takes to radiate away its total initial mass Mq. Easily, we get the expressions 
for the Tev corresponding to (|66|), and (|67|), concerning the values of the constant parameter C in (|7^) and j?^), 
respectively. Thus, they can be written as 



and (75) 




Te^n = 7 \ / . r \ • (76) 



It is clear that given the initial values for the mass and radius (therefore, the parameter C = -j^ ), Tev for both 
evolutions is fully determined. It is also worth mentioning the particular astrophysical scenario \h&t the evolution 
of these matter configurations could simulate. Both of them (either (|66|), or (|67|)) start at a particular value of the 
parameter C, completely radiate away their initial mass and disappear leaving no remnant. In this sense this scenario 
resembles Type I Supernova but, in the present case, it begins from a neutron star progenitor (characteristic values 
for the radius, R% and the mass Mq) and gravitation plays an important role. In the standard Type I Supernova 
picture there is a white dwarf progenitor and gravitation is not longer the main cause of the emission processes. 

In Figure 1 the total evolution times Tevi and Tevu are sketched as a function of the initial mass. In both cases 
we have assumed 10 Kms. as the value for the initial radius. This value can be considered as typical for neutron 
stars 10,^ ■ It can be appreciated from this figure that as the initial mass Mq rises, the total time, Tevj, increases 

asymptotically up to the value Mq — . The other branch Tevjj decreases coming from this same asymptotic 
maximum for Mq. These two branches of the evolution time may represent different astrophysical phenomena. In 

first case, Mq < , it is apparent that the more massive the configuration is, the more time it takes to radiate 

the initial mass. In the second branch, Mq > more massive spheres lead to more violent collapses. The value 

of Mq = emerges as a critical point of stability for a configuration with a NLES, having a CKV and matched to 
a Vaidya exterior solution. Any perturbation accreting over (or expelling mass from) this type of relativistic sphere 



will initiate processes which radiate away a significant amount of energy. 

In the next section we shall consider the hydrodynamic consequences of configuration satisfying a NLES for both 
cases: C =const. and C = C(T). 



V. COLLAPSING RADIATING SPHERES 



Finally, to determine if anisotropic fluid spheres with a NLES can conform a feasible matter candidate to stellar 
models, we shall study collapsing scenarios for radiating configuration of these type of materials. In what follows, 
radiating collapsing models will be matched to the Vaidya exterior metric and energy density will always be positive 
and larger than pressure everywhere within the fluid distribution. Finally, the fluid velocity, as measured by the 
locally Minkowskian observer, will always be less than one. 

Besides these "regularity" conditions, the modeling emerges from a heuristic assumption relating density, pressure 
and radial matter velocity. This ansatz, guided by solid physical principles, allows the generation of radiating solutions 
from a known static "seed" solutions and reduces the problem of solving Einstein Equations to a numerical integration 
of a system of ordinary differential equations for quantities evaluated at the surfaces (shocks and/or boundaries). The 
rationale behind this ansatz can be grasped in terms of the characteristic times for different processes involved in a 
collapse scenario. If the hydrostatic time scale, Thydr, which is of the order ~ 1/^/Gp (where G is the gravitational 
constant and p denotes the mean density), is much smaller than the Kelvin- Helmholtz time scale {Tkh), then in a 
first approximation, the inertial terms in the equation of motion can be ignored p6[. 
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A. The metric and the matter 



The method we use to build these models is a general strategy presented several years ago by L. Herrera, J. Jimenez 
and G. Ruggeri (HJR). Only a very brief description of this method is given here. We refer the reader to |17| and 
[ p^ for details. In this method a nonstatic spherically symmetric distribution of matter is assumed. The metric 
representing this distribution, in Bondi radiation coordinates p9| , takes the form 



z'^Pi^^^)h{u,r) du^ + 2e2'^("''') dudr - r^{dd^ + sh 



with h=l 



2m [u, r) 



(77) 



Here u — x'^ is a. time like coordinate, r — is. the null coordinate and Q ~ x^ and tp = x^ arc the usual angle 
coordinates. The u -coordinate is the retarded time in flat space-time and, therefore, u -constant surfaces are null 
cones open to the future. The function m (w, r) is the generalization, inside of the distribution, of the "mass aspect" 
defined by Bondi | [Tg| ] which in the static limit coincides with the Schwarzschild mass. In addition we consider an 
anisotropic fluid sphere composed by a material medium plus radiation. For a observer with a radial velocity —uj with 
respect to the fluid, the energy-momentum tensor can be written as |18|,p(| 



= {P + Pi 



- P^^t^y + {Pr - PL)n^n^ 



(78) 



with 



(1 



hi 



(1 



hi 



-51 



(79) 



and where the physical variables, as measured by a local minkowskian observer, are represented by the energy 
density, p; the radial and tangential pressure, Pr, P±, respectively; the radiation energy flux density and fluid 
radial velocity, ~uj. The co- moving minkowskian observer coincides with the Lagrangean frame (the proper frame) 
where the interaction between radiation and matter can be easily handled Q . Notice that the radiation is treated in 
the diffusive regime, therefore it is considered to have a mean free path much smaller than the characteristic length 
of the system. Within this regime, radiation is locally isotropic and we have 



T = q 



(80) 



Where, P, describes the fluid radial pressure and the radiation contribution to the energy density, energy flux density 
and radial pressure are represented hy pn, J- and V, respectively. It is clear that in this radiation limit the total radial 
pressure, Pr, encompasses the hydrodynamic and the radiation contribution to the pressure (see for details). 



B. The Einstein Field Equations and the Matching 



Inside the matter distribution the Einstein field equations can be written as 

p + 2ujq + UJ^Pr 1 



47rr2 \ he^P 
■ UjPr — (1 — uj)q 1 



1 -f 47rr^ 
{p + Pr - 2q)[l - UJ 



Pi 



47re2/3 



1 -|- a; 2'Kr 
3/3'(l - 2m') - m" 
87rr 



0, 
0, 
0, 
= 0. 



and 



(81) 
(82) 
(83) 
(84) 



As was stressed by Bondi |l^ it is possible to solve algebraically the physical variables of equations (81) - (|8J) if we 
have both m (u, r) and (3 {u, r) (and their derivatives) plus another equation that relates the tangential and radial 
pressure. This method will be clarified by an example in the next section. 

From the field equations (|8|) and (pSh it is clear that they can be (formally) integrated yielding 
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m — 4:n 



■ toPr — (1 — (^)q 

l+U! 



r^dr = An I r^pdr , 




and 



' l+u h 



where 



p — ujPr — (1 — Lu)q 

l+U! 



and 



P 



p + P] -dr, 



-ujp + Pr - (1 - uj)q 

1+UJ 



(85) 
(86) 

(87) 



These auxihary functions p and P must coincide with the energy density and the radial pressure in the static hmit. 

Matching the Vaidya metric to the Bondi metric (77) at the surface of the fluid distribution, r — a, implies /3a = 
with the continuity of the mass function m (u, r), i.e. the continuity of the first fundamental form. The continuity of 
the second fundamental form leads to (again, see [H for details) 

Pa ^ Pa ■ 

The subscript a indicates that the corresponding quantity is evaluated at the boundary surface. 



C. Surface equations 

As we have mentioned, the crucial point of the HJR method is the assumption that the radial dependence of p 
and P can be borrowed from a "seed" static solution and, therefore the metric functions m (u, r) and /? (m, r) can be 
determined from eq. ( ^5| ) and (|8^), up to some functions of the time-like coordinate u, which in turn are computed 
by integrating a system of ordinary differential equation coming from the boundary conditions. This assumption 
represents a correction to the first approximation in the radial velocity (i.e., quasi-stationary approximation) and is 
expected to yield good results whenever Tkh 3> Thydr- Fortunately enough, this condition is fulfilled for almost all 
kind of stellar objects. For example, in the case of the Sun we get Tkh ^ 107 years, whereas Thydr 27 minutes. 
Also, the Kelvin- Helmholtz phase of the birth of a neutron star last for about tens of seconds, whereas for a neutron 
star of one solar mass and a ten kilometer radius, we obtain Thydr ~ 10~^sec. 

It is clear that with the metric functions completely determined, the physical variables (p, Pr, P± , and uj) are 
algebraically solved from the field equations (^ij) through (^) . 

Scaling the radius a, the total mass ma and the timclikc coordinate u by the total initial, mass ma (0), i.e. 



A{u) 



a{u) 
nia (0) ' 



M{u) 



ma (0)' 



and defining 



, , 2M(u) 
F{u) = l ttV, and 



and u 



ma (0) ' 



1 



(89) 



(90) 



A{u) ' ' ' 1 - cjq 

The first surface equation comes from the definition of the velocity in radiation coordinates, as 

A^F{n~l). (91) 
The second surface equation arises from the luminosity evaluated at the surface which can be written as 

L = -M = inA^qa (2^2 - 1) F . (92) 



Using equation (91) and the definitions (|9C|), we can re-state equation ( |92|) as 

2L + F{1-F) (n-l) 



F = 



A 



(93) 



Finally, some straightforward manipulations coming from the field equations (p2|), (p2[), ( |83| ) and (§J), lead to 




R--(Pr-P] =0, 



(94) 
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where 



R = P' 



p + P 



AnrP 



iP± - Pr) 



(95) 



Notice that (joj) corresponds to a generaUzation of the Tolman-Oppenheimer-Volkov equation for any dynamic radia- 
tive situation. Evahiating it on the boundary surface r = a- and using the definitions (pO|), the third surface equation 
can be obtained, i.e. 



F 
F 



Pa 
Pa 



2Fn Pr, 



Pa 



A 



Pa 



G = 0, 



where 



At: A (3^2 - 1) , 3 + F 
2A 



-Pa 



Fn^^'-^iP. 

Pa Apa 



Pr 



(96) 



(97) 



Equations (91), ( |93| ) and ( |9q ) constitute the System of Surface Equations (SSE). This system can be integrated for 
any given radial dependence of the effective variables p and P. In the context HJR with anisotropy a general equation 
has to be provided in order to relate the tangential and the radial pressure |23Jl^,^ . 



P± (u, r) - Pr {u, r) 



rP' 



p + P 



(98) 



It is clear that if the anisotropic equation of state ( |98| ) is given, we have three differential equations to get four 
unknown functions, i.e. the exterior radius A{u), the gravitational potential at the surface, F{u), the radial velocity, 
^(u), and the luminosity profile, L{u). Since the only observable quantity entering a "real" gravitational collapse is 
the luminosity, it seems reasonable to provide such a profile as an input for our modelling. 



D. The HJR formalism and NLES 



At this point we shall verify that, in the HJR formalism, models having 
within the distribution. Taking derivatives of (|^) with respect to r, we get 



~l3'e-^f^C{u), 



and using 



) and (86) with (M) we obtain 



m — 2'Kr' 



satisfy a nonlocal equation (|| or H) 

(99) 

(100) 



On the other hand, taking derivatives of (IOC) with respect to r and comparing it with ( p5[ ) it follows that 

p-W + r{p' - P'^ = Q. (101) 
As expected, this expression emerges from the condition ^ and in terms of the field equations it reads 

G« + 3G; + r[(G:;i)' + (G;:)'] =0. (102) 



It should be noticed that from equation (101) two expressions for NLES can be obtained. They can be written as 

(103) 
(104) 



2 f H(u) 
p{r,u) = P{r,u) + - / P(f, u)df+ ^ ' 



or 



P{r,u) — p{r,u) g / r p{r,u)dr 

Jo 



r 

B{u) 



In the present case, equations (103) and (104) are (||) and (^, respectively, expressed in terms of the effective variables 
defined in (p^). From this new aspect of the NLES, it is clear that the contribution to the pressure and the energy 
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density coming from the radiation energy flux, as well as the effects of the fluid velocity, are taken into account in the 
nonlocal character of the equation of state. Moreover, equation (^), in terms of the effective variables, becomes 

P(r, u) = Ipir, u) - 2 (lp{r, u) - I (p(r, u))]+^. (105) 



Clearly, it is the dynamic generalization of equation (|1^), with the corresponding interpretation of P{r, u) — ^p{r, u), 
as the "effective equation of state" which may be affected by the standard deviation 

c^ip(r,«) = (^^Pir,u) - ^ {p{r,u))^ , (106) 

coming from the contribution of the nonlocal term (/5(r, u)). 

There is an other important consequence arising from (g) in terms of the surface variables in the HJR formalism. 
First, comparing the condition (^ with the definition for F {u), it follows 

ha = l-2^ = F{u)^C{u). (107) 



Notice that C(m) is related to the gravitational potential at the surface. Next, evaluating (IOC) at the surface and 
using ( p8|) we get 

rUa = 2T:a^ [1 + LUa) Pa ■ (108) 

This equation can be written as a function of the dimensionless variables (B9|) as 



l-C{u)= {2n{u)-l)pa{u) , (109) 

\l (u) 

therefore, the "velocity" of the surface (u) can be obtained 

^ AnAjufpAu) 

8TrA{uf pa{u) + C{u)~l 

It is clear that it relates the velocity and the effective density at the surface. 

Now we are going to provide several "seed" solutions to explore reasonable collapsing scenarios compatible with a 
NLES in the HJR formahsm. 



E. Collapsing Scenarios 



The first example considered as a starting static equation of state is the Schwarzschild solution. This model 
represents a generalization of a radiating incompressible anisotropic fluid of homogeneous density. The effective 
density can be written as 



Psh 



1 

8^ 



/(") 



(111) 



With equation ( p5[ ) evaluated at the surface we find f (u) and using (IIC) we obtain that the fluid velocity at the 
surface is constant and takes the value Wa = — |c . 

The second case of study is the anisotropic Tolman-VI-like model. Again, in the static limit this solution is not 
deprived of a physical meaning. The static Tolman VI solution approaches that of a highly relativistic Fermi Gas 
and, therefore, one with the corresponding adiabatic exponent of 4/3. In this case we have 



1 



PTvi 



Sirr 



(112) 



Again equation (85), can be evaluated at the surface, with the restriction (110), leading to a constant velocity at the 
surface: uja = c. 
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Although the above equations of state have been successfully worked out within HJR formalism, in our study with 
a NLES they lead to very restrictive matter configurations where the velocity of the surface is held constant during 
the collapse. It is worth mentioning that the above results for the Schwarzschild-like and Tolman- Vl-like models are 
valid for the general case C — G{u). 

Next, we use as a "seed" another static equation of state which represents a richer description for ultracompact 
nuclear matter. We shall explore the consequences of adopting a static solution proposed by M. K. Gokhroo and 
A. L. Mehra in 1994 This solution corresponds to an anisotropic fluid with variable density. It has been used 
in the context of the HJR method to study several interesting properties of an anisotropic flui d |^0[| describing the 
Kelvin-Helmoltz phase in the birth of a neutron star |2^,^. It leads, under some circumstances |20 , to densities and 
pressures given rise to an equation of state similar to the Bethe-Borner-Sato newtonian equation for nuclear matter 

The energy density and radial pressure for this static space-time are assumed to be 



p{r) 



P{r) 



Pr 1 



, (0<A:<1) 



2m (r) 



and 
with 



n> I: 



(113) 
(114) 



where pc, Pc and a are the energy density, radial pressure and the surface of the sphere, respectively. Finally, k is an 
arbitrary constant which may take values between zero and one and a constant parameter n > 1. 



F. The modelling performed 



In applying the HJR method to the Gokhroo and Mehra model we assume that the effective variable p has 
the same r dependence as p (r), i.e. 



p{u,r) = Pc (u) 



1 - K{u)- 



„21 



, (0<if(u)<l) 



(115) 



In the above equation, K iu) is a function that may take values between zero and one only |20| , and the energy density 
at the center of the configuration, pc (u) , coincides with the pc in the static case. 

It is necessary to relate the unknown functions K (u) and pc (u) to the surface functions A{u), fl(u) , F(u) and 
L{u). Thus, we evaluate the expression (p3) at the surface of the configuration: 



m„ = M - 



— (5 - 3K{u)) = - (1 - C (w)) . 



Therefore, from the above equation it is possible to find '■ 

, _ 15(l-C(u)) 



87rA2 {5-3K{u)) ' 



and the expression for K (u) is obtained from (110), namely 

, , 5 417 - 3 5 1 + 3 Wa 
K [u) = - ■ 



3 817-5 33 + 5 Wa 



Thus, the effective energy density is written as 



p{r,u) 



[l-C(M)]r2 
167rA4 (217 - 1) 



42 

5 (3 - 417) + 3^ (817 - 5) 



and from ( |104| ) the effective pressure, can be computed as 



P{r,u) = 



[1 - C (u)] r2 



A^ 



3 (3 - 417) + — (817 - 5) 



(116) 



(117) 



(118) 



(119) 



167rA4 (217 - 1) 

Notice that the function B{u) should vanished in order to fulfil the matching condition (|88|), i.e. B{u) = 0, Vu. 



(120) 
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Finally, by virtue of 

m {u, r) 



and (|8^) the metric functions read 
r3(l-C(u)) 



2'm{r, u) 



and 



1 



(121) 
(122) 



Next section is devoted to studied the p hysic al co ntrib ution of the NLES to several collapsing model s em erging from 
the above effective variables (equations (119) and (|120[)) and from the corresponding metric functions (121) and (122). 



G. The Systems of Surface Equations 

At this point two family of models will be considered. The first family emerges with C — const, which, by equation 
(107), clearly represents a configuration with a constant gravitational potential at the surface, i.e. 



1 



F ~ C ~ const. 



(123) 



The second family of models is the most general one. It is assumed that C is a function of the time like coordinate, 
i.e. C = C{u). 

In order to understand the contribution of the NLES to the hydrodynamic scenario during the collapse, we shall 
compare our models with NLES to the most similar standard Gokhroo-Mehra-Martfnez model ||2^] . 
In the first case, C = const, the system of surface equations is 



L = -ic(17-l)(l-C) 



n = 



2n- 1 



L-i(l-C) 



and 

n- 1 



2A 



[1-11C-4C (8rj-9) n] 



(124) 
(125) 

(126) 



This system of nonlinear ordinary differential equations can be solved (numerically) for a given set of initial values 
for A{u), F (u) and (u). 

As we have stressed above, the reference model will be the Gokhroo-Mehra-Martfnez model with a constant 
gravitational potential at the surface of the configuration. In this model the only di ffere nt equation is the third surface 
equation. Thus, the system of surface equations is conformed by equations (124), (125) and by the equation (4.28) in 
reference |2^] assuming F — 0. Notice that, in this case, the luminosity profile, L (u), emerges as a consequence of the 
condition that the gravitational potential at the surface of the configuration be constant. 

In the case C = C(m), the system of surface equations is given by 



C(m 



A = c(u) [f^ - 1] , 

2i + C(w) [1 - C(m)] [17 - 1] 



A 



L-A{1- C{u)) 



and 

ACju) 
" 2C(u) 



n - 1 

2A 



[1 - 11 c{u) - 4C(it) {sn - 9) n] 



(127) 
(128) 

(129) 



The set of initial condition and parameters for three models are 



Af(0) = 


I.OA/q 


Pc(0) = 


2.2 X 10^* 


g cm 


A{Q) = 


10 


Pa(0) = 


9.8 X 1013 


g cm~^ 


n{o) = 


1 - e 


a(0) = 


1477 m 




K{0) = 


5/9 


ZaiO) - 


0.118 




C(0) = 


4/5 


UJa (0) = 


-1.0 X 10" 


-10^ 



Concerning the above set, it should be pointed out that: 
• All of them correspond to typical values for young neutron stars. 
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• Because of (^0|), £7 = 1 corresponds to = . From the above systems of surface equations, it is straightforward 
to verify that this initial velocity lead to static models. 

• In order to obtain collapsing configurations the initial surface velocity has been perturbed with an e fvlO^^'^. 

• In our simulations, we have imposed that the energy conditions for imperfect fluid be satisfied. In addition, the 
restrictions — 1 < w < 1 and r > 2m {u,r) at any shell within the matter configuration are also fulfilled. 

• For C = C(u) : 

— we do not need any perturbation in the velocity to start the collapse. 

— The luminosity profile, L (u), has been provided as a Gaussian pulse centered aX u = Up 

with A the width of the pulse and AMrad the total mass lost in the process. This family of models has 
been simulated using 

AMrad = 10"^ M (0), A = 5, Up = 100 ms. (131) 



H. The Evolution of the physical variables 



The evolution of the physical variables are presented as functions of the standard Schwarzschild time. The rela- 
tionship between the usual Schwarzschild coordinates, (T, _R, 6,$), and the previously mentioned Bondi's radiation 
coordinates can be expressed as 



u = T- ^737 7r-( AT dr, e^e, r = R and = $ . (132) 

J e^P[r — 2m[u, r)) 

Figure 2 displays the evolution of the boundary of the surface for the Gokhroo-Mehra model with a NLES (plate (a)) 
and with constant gravitational potential at the surface (plate (b)). None of them tend to any asymptotic equilibrium 
configuration as was reported earlier . It is apparent from these figures that our models with a NLES lead to a 
softer configurations than those without this particular equation of state. Our matter configuration starts to collapse 
earlie r and evolves deeper than the reference Gokhroo-Mehra-Marti'nez model with F = 0. In addition, because of 



(123) M (X A , the models with a NLES radiates more energy than those without it. 

The profiles of the physical variables (p, Pr, q , and uj/c ) for three different times, are sketched in Figures 3 and 
4 (plates (a) through (d)). Figure 3 displays the profiles for the Gokhroo-Mehra model with a NLES and Figure 
4 Gokhroo-Mehra-Martinez model with constant gravitational potential at the surface. It can be appreciated from 
these figures that the changes in time of the density and pressure (Figure 3, plates (a) and (b)) of the configuration 
with the NLES are more significant than those observed in Figure 4, (plates (a) and (b)). In fact, from Figure 3 (b) 
it is clear that the pressure rises at the inner core while it diminishes at the outer mantle. This can be understood if 
we recall that the total radial pressure Pr collects the contribution of both the hydrodynamic and radiation pressures 
and that it is at the outer layers where the radiation pressure becomes more significant in time. Now, concerning 
plates (b) and (c), we can see that around the zone of r = 5.5 Kms., the total pressure remains constant and coincides 
with the maximum of the energy flux density. It is clear that the energy flux density in models with a NLES (Figure 
3, plate (c)) is more intense than that observed in the Gokhroo-Mehra-Martfnez reference model (Figure 4, plate (c)). 
Finally notice that in Figures 4(c) and 4(d), those shells that have a strong energy flux density also have the higher 
speed. This situation contrasts with the difference that can be observed in the Figure 3, (plates (c) and (d)) between 
the maximum of the energy flux density (around r — 5.5 Kms.) and the maximum in the velocity at the r = 9.2 
Kms. This difference of about 4 Kms. remains constant in time during the collapse. This effect also appears in the 
NLES with variable gravitational potential at the surface (Figure 5, plates (c) through (d)), and is not present in 
the corresponding Gokhroo-Mehra-Martfnez reported in ]2C| ]. It is worth mentioning that, because of our assumption 
of the diffusion regime approximation for the radiation flux, we expect that hydrodynamics and radiation should be 
closely related and this is clear in the Gokhroo-Mehra-Martfnez models. Therefore, it seems that the shift between 
the radiation maxima and the velocities could be thought of as the main effect emerging from the integral contribution 
presented in the nonlocal equations of state (^) and (||). 
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Finally, in plates (e) and (f) in Figure 4 particular hydrodynamic discontinuities, that should be mentioned, are 
displayed. These plates zoom in time of the corresponding plates (c) and (d) of the same figure. As it can be 
appreciated from these plates, as the energy flux rises there is a wave front propagating outward. The discontinuity 
is present not in the physical variable (either the energy flux or the matter velocity) but in its derivative. This wave 
front, driven by the radiation pressure, can be associated to a hypersurface traveling with the velocity of sound [ p6[ . 
It is clear that the shell where the discontinuity of the derivative of the energy flux takes place coincides with the 
corresponding surface where the discontinuity in the derivative of the velocity is identified and this coincidence is 
consistent with the radiation regime that we have assumed. 



We have shown that fiuids satisfying nonlocal equations of state (either (|5|) or (|6|)) could represent feasible bounded 
matter configurations and some of them could have significant geometric properties. 
Some answers to the questions that motivate the present work have been explored: 



• Assuming (Efl) and C as a constant parameter in (0), we have found that spherically symmetric space-times 



• The anisotropic imperfect fluid, emerging from the metric (|3^), fulfills the corresponding energy conditions. 

• Bounded gravitational sources can materialize from this type of fluids and meaningful hydrodynamic scenarios 
can be obtained. 

• Although the simulations for collapsing configuration where carried out for a very short period of time, it seems 
that the radiating spheres we have considered could represent some initial phases in the evolution of compact 



More investigations on this curious equation of state are needed to extend our understanding of the possible collective 
effects that are clear in cither (|^) or (^. But, as the main conclusion of the present work it can be asserted that 
Nonlocal Equations of State could represent satisfactory candidates to describe the behavior of matter at supranuclear 
densities in General Relativity. 
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VI. CONCLUDING REMARKS 




objects. 
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Figure Captions 



Figure 1 The total evolution times Tevi and Tevn are sketched as function of the initial mass. In both cases we 
have assumed 10 Kms. as the value for the initial radius of the collapsing configurations. 

Figure 2 Evolution of the boundary radius of the matter configuration satisfying NLES, (plate (a)) and for the 
standard Gokhroo and Mehra model with constant gravitational potential at the surface constant gravitational 
potential at the surface (plate (b)). 

Figure 3 Evolution of matter variables for the Gokhroo and Mehra model with a NLES. Plates (a) through (d) 
represent the profiles for the density, pressure, energy flux density and mass velocity, respectively at three 
different times. 

Figure 4 Evolution of matter variables for the standard Gokhroo and Mehra model with constant gravitational 
potential at the surface. Plates (a) through (d) represent the proflles for the density, pressure, energy flux 
density and mass velocity, respectively at three different times. Plates (e) and (f) are a zoom in time to the 
corresponding plates (c) and (d), respectively. 

Figure 5 Evolution of matter variables satisfying NLES with F{u) = 1 — ^^^^ = C{u). Plates (a) through (d) 
represent the profiles for the density, pressure, energy flux density and mass velocity, respectively at three 
different times. 
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